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Abstract 

We compute the leading-log QED, the next-to-leading-log QED-QCD, and the 
electroweak corrections to the charm quark contribution relevant for the rare de- 
cay K+ — » Tt + vv. The corresponding parameter P C (X) is increased by up to 2% 
with respect to the pure QCD estimate to P C (X) = 0.372 ± 0.015 for m c (m c ) = 
(1.286±0.013) GeV, a s (M z ) = 0.1176±0.0020 and \V US \ = 0.2255. For the branching 
ratio we find B(K + — > -n + vv) = (8.5 ± 0.7) x 10 -11 , where the quoted uncertainty is 
dominated by the CKM elements. 

1 Introduction 

The rare decay K + — > 7i + uu is both theoretically very clean and highly sensitive to short- 
distance physics and thus plays an outstanding role among flavour-changing neutral current 
processes both in the standard model (SM) and its extensions [1-3]. Together with the 
process Kl — > ir vv it provides a critical test for the Cabibbo-Kobayashi-Maskawa (CKM) 
mechanism of CP violation, while it probes operators generated by new physics at energy 
scales of several TeV [4]. 

In the SM the decay K + — > -n+vv proceeds through Z-penguin and electroweak box 
diagrams of 0(G 2 F ) which exhibit a power-like GIM mechanism. This implies that non- 
perturbative effects are severely suppressed and, related to this, that the low-energy effec- 
tive Hamiltonian [5,6] 

H cS = a ( X cX l (x c ) + X t X{x t )) {s Ll ,d L ){v lLl »v lL ) (1) 

l=e,fi,T 

involves to an excellent approximation only a single effective operator. Here Gp is the Fermi 
constant, a the electromagnetic coupling and 6w the weak mixing angle. The sum is over 
all lepton flavours, Aj = V* s Vid comprises the CKM factors and fa represents left-handed 
fermion fields. 

The function X(x t ), where x t = Tn^(fit)/M^ and mf(fi t ) is the top quark MS mass, 
describes the matching contributions of internal top quarks to the operator in Eq. ((Tj), where 
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Figure 1: Examples of leading-order diagrams contributing to the decay K 
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7r + z^ in the SM. 



the matching is carried out at the scale fit = 0(m t ). Sample diagrams are shown in Fig. [U 
The energy scales involved are of the order of the electroweak scale or higher, while both the 
QCD and QED anomalous dimensions of the corresponding operator vanish. Hence X(x t ) 
can be calculated within fixed-order perturbation theory. The relevant Z-penguin and 
electroweak box diagrams are known through next-to-leading order (NLO) in QCD [6-9]. 
The inclusion of these 0(a s ) corrections allowed to reduce the ±6% uncertainty related to 
the top quark matching scale fit = 0(m t ) present in the leading order (LO) formula down 
to ±1%. The leading term in the large top quark mass expansion of the electroweak two- 
loop corrections typically amounts to a per mil correction for the branching ratio if the MS 
definition of a and sin 2 Q w is used, while the uncertainty related to unknown sub-leading 
electroweak contributions is conservatively estimated to be ±2% [10]. 

The function X e (x c ), relevant only for K + — > 7i + uu, depends on the charm quark MS 
mass through the parameter x c , conventionally defined as 

As now both high-energy and low-energy scales, namely [iw — 0{M W ) and fi c = 0{m c ) 
are involved, a complete renormalisation group analysis of X e (x c ) is required. In this 
manner, large logarithms ln(/i 2 / /i 2 ^) are summed to all orders in a s . At LO such an 
analysis has been performed in [11]. The large scale uncertainty due to /i c of ±26% in this 
result was reduced by a NLO [5,6] and a subsequent NNLO calculation [12-14] to ±2.5%. 
While the QCD part of the calculation has reached a high level of sophistication no QED 
or electroweak corrections have been included so far. We close this gap by calculating 
the LO and NLO logarithmic QED corrections as well as fixing the scheme of the input 
parameters in sin 2 9w and a by an electroweak matching calculation. The latter point 
can be exemplified by noting that the charm quark contribution is mediated by a double 
insertion of two dimension-six operators. This results in a contribution of 0(G 2 F ) - the 
second power of Gp resides in x c - plus electroweak corrections. Yet the leading result 
of Eq. ([T]) can only approximate the electroweak corrections for a specific choice of the 
renormalisation scheme for the prefactor of the charm quark contribution, expressed as 
a/ sin 2 6w- While it is expected that using MS parameters renormalised at the electroweak 
scale would approximate the electroweak corrections best [15] only an explicit calculation 
can provide a definite result. In this work we normalise all dimension-six operators to Gp. 
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Thus, we replace the parameter x c in Eq. with the unfamiliar definition 



V2 



sin 2 8 W 
ira 



G F m 2 c (fi c 



(3) 



which only at tree level equals the familiar ratio m 2 (/i c )/M^/. 

The hadronic matrix element of the low-energy effective Hamiltonian can be extracted 
from the well-measured Ki 3 decays, including isospin breaking and long-distance QED ra- 
diative corrections [16-18]. After summation over the three neutrino flavours the resulting 
branching ratio for K + — > n + vv can be written a£] [5,6, 19] 
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ImA/ 
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ReA c 



(P C (X) + 5P C , U ) 



ReXi 



X(x t ) 



(4) 



The parameter 
Pc(X) = 



1 

A 1 



-X e 
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(5) 



describes the short-distance contribution of the charm quark, where A = \V US \. The charm 
quark contribution of dimension-eight operators at the charm quark scale /i c combined 
with long distance contributions were calculated in Ref. [19] to be 



SPr 



0.04 ±0.02 



(6) 



The quoted error on this value can in principle be reduced with the help of lattice QCD [20]. 

The remaining long distance corrections are factored out into the following two param- 
eters: k + contains higher-order electroweak corrections to the low energy matrix elements, 
and A EM denotes long distance QED corrections. A detailed analysis of these contributions 
to NLO and partially NNLO in chiral perturbation theory has been performed by Mescia 
and Smith in [17], who found the numerical values k + = (0.5173±0.0025) x 10~ 10 (A/0.225) 8 
and A EM = -0.3%. 



2 Electroweak Corrections in the Charm Sector 

The charm quark contribution involves several different scales and the corresponding large 
logarithms have to be summed using renormalisation group improved perturbation theory. 
Keeping terms to 0(a s ) and 0(a/a s ) the expansion of the parameter P C (X) reads 

P C (X) = ^P^\X) + P^\X) + ^PP\X) + 4^P!:^X) + ^P^ (7) 
a s (ii c ) 4vr aiWc) ots(Hc) 

have omitted a term which arises from the implicit sum over lepton flavours in P c because it 
amounts to only 0.2% of the branching fraction. 
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The LO term P C (0) (X), the NLO term P C (1) (X), and the NNLO term pP(X) have been 
calculated in [11], in [5,6], and in [14] respectively. The main goal of the this paper is to 
present the electroweak corrections P c (X) and Pc £S \x). 

The calculation is performed in two steps. First, at the scale nw ~ M w the SM is 
matched to an effective theory where the top quark, the W boson, and the Z boson are 
integrated out, but the charm quark is still a dynamical degree of freedom. Second, at the 
scale /i c ~ m c the charm quark is integrated out and the effective Hamiltonian in Eq. (CQ) 
is obtained. 

After integrating out the particles at the electroweak scale the effective Hamiltonian 
containing the dimension-six operators takes the following form: 

4G f 



q=u,c 



+\ c J2cAv)(Qj - Qj) + \c A (fi)Q A + \c v ^)Q, 

Here we kept only operators relevant for the decay K + — > 7i + uu. These are the semi- 
leptonic operators 

Q3q = ^2 (s L -f^q L )(ui L ^l L ) and Q 4q = ^ ^Ll^d L ){l L ^v lL ) , (9) 

the current-current four-quark operators 

Ql = \ {{stl»QlWh<) ± {sllAMl'dl)) , (10) 
where a, (3 are colour indices, and the operators 
Qi = EE (-^ 3 )(<?757^)(^7^l) , 

(11) 

Qv = Y,zZ (I 3 q -^Qq^ 2 0w)(qi,q)(^Ll^iL), 

q l=e,/j,,T 

which describe the quark-neutrino interaction. We follow Ref. [14] in the definition of the 
evanescent operators. All evanescent operators relevant first at the order considered in this 
work are defined as 

E S = E (^T^T^X^VV 3 ^) - (16 - 4e)Q 3g , (12) 

l=e,fi,r 

i.e. the evanescent operator needed for the QED renormalisation of Qsq, or in an analogous 
way. 

These operators mix via double insertions into the operator given in Eq. (CQ). Tradi- 
tionally one distinguishes the box contribution which comprises double insertions of the 
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Figure 2: Leading-order diagrams for the mixing of various dimension-six operators into Q u (see 
text for details). 

semileptonic operators Q 3q and Q ig (see Fig. El left side) and the penguin contribution 
which comprises double insertions of the current-current-type operators Q± and the oper- 
ators Qa and Qz (Fig. El right side). The relevant dimension-eight part of the effective 
Hamiltonian can then be written as 

^charm = ^G%X e C^ (jl) + (? F \P* (jlj) Q v , (13) 

where the operator Q v is defined as 

2 

Q» = E ^liM^ilI^il) , (14) 

while and denote the box and penguin contribution, respectively. 

The renormalisation group analysis proceeds in several steps. The initial conditions 
for the renormalisation group equations (RGE), which govern the running of the Wilson 
coefficients, are calculated in Sec. 12.11 The anomalous dimensions are computed in Sec. 12.21 
After integrating out the bottom and the charm quark, the theory is matched onto the low 
energy effective Hamiltonian of Eq. (pQ). The relevant results are collected in Sec. 12.31 In 
Sec. 12.41 the pieces are put together to give the final result for P C (X). 

We have computed all Feynman diagrams in this paper using FORM [21] routines and 
independently using Mathematica. All the QCD corrections relevant to a NNLO analysis 
of P C {X) are given in [14] and references therein. 

2.1 Initial Conditions 

The Wilson coefficients are found by matching the one light particle irreducible Green's 
functions in the full and the effective theory at the electroweak scale \x\^ ~ M^. We use 
the MS scheme for both theories and remark that a finite field redefinition for the light 
particles ensures the correct normalisation of the kinetic term in the effective theory. In 
the box sector only Cw and in the penguin sector only C± and Ca/v receive electroweak 
corrections at the order considered here (see Fig. [3]). We expand the Wilson coefficients in 
powers of the coupling constants 

C(ji) = C^(fx) + ^Mc^(fi) + -C^(fi) + -^c (es V) (is) 

47T a s 47T 
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Figure 3: Feynman diagrams contributing to the NLO matching for Ca- 



and use a similar expansion for any quantity in the following, unless explicitly stated 
otherwise. 

We normalise the Wilson coefficients Cw, C±, and Ca/v to the muon decay constant 
Gp [22] . In this way most of the radiative corrections cancel, including all terms dependent 
on m t and Mh in the case of Cy/ and C±. All our matching calculations have been 
performed in the generalised gauge for the photon field and in the case of Ca also for 
the W and Z fields as a check of our results. 

At the one-loop level a neutrino-photon Green's function is generated which contributes 
to Cy via the equations of motion. Yet Qy does not mix into Q u and the Wilson coefficient 
Cy is not needed. 

For the relevant electroweak corrections we find 



C±{nw) — 1 > 

4 e W) = o, 

in agreement with Ref. [23,24], 



(16) 



Cw(^w) — , 
<#W) = -£-21»* 



and 



C A °\nw) — 1 , 

(es) _ 3m f 2 Hag, - 6 3 - c|Mf Mg, 
, 3M| Mf 



(17) 



4(M^-M^)(M^-cX) M; 
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Figure 4: Sample two-loop diagrams contributing to the self mixing of Qs c - Wavy lines denote 
photons, curly lines denote gluons. 

2.2 Anomalous Dimensions and RGE 

The mixing of dimension- six into dimension-eight operators through double insertions leads 
in general to inhomogeneous RGE [25]. In the box sector they are given by 

t*j-C?(n) = lv C*(n) + A^C w {n)C w {n) , (19) 

li—C w {n) = j w C w {fJ,) , (20) 
an 

where is the anomalous dimension of Qs q , ^ v encodes the running of Q u , which stems 
solely from the running mass and coupling constant which in our definition multiply the 
Q u operator, and 7^ is the anomalous dimension tensor of the mixing of the operators Q% q 
and Q^ q into Q v . 

7^ is given in terms of the QCD /3-function and the anomalous dimension of the charm 
quark mass by 

7^ = 2(7^-^). (21) 
The explicit values are 

(o) =8 ~{*) = l 7 («) = ^ (22) 

Im u > Im g ' Im g ) V / 

A)=ll-?/, & = 0, ft. = -jj(/„ + j), (23) 

where / u and denote the number of up- and down-type quark flavours, and / = f u + /<j. 

The remaining anomalous dimensions can be calculated from the pole parts of one- and 
two-loop diagrams, some of which are shown in Figs. H] and 0, using standard methods 
[25-27]. We find the following values: 

7f 0) = "8, 7 f e) =0, 7 f es) = -^, (24) 
7$ = 0, 7$ = -4, 7^ = 4- (25) 
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Figure 5: Sample diagrams for the NLO mixing of Q% q and Qt±q into Q t 
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B(0) 



is known for a long time (see [5] and references therein), and 7^ and 7jy Sj have 
already been calculated in [22]. 

In order to solve the RGE we perform a trick [5,28], so that we can use the RGE for 
single insertions also in our case. To this end, we rewrite Eq. ( 1201) as 



Then we can combine both equations (fT9i) and ( 1201) into a linear equation 
d 



(26) 

(27) 
(28) 

(29) 
(30) 



where 



The RGE for the penguin sector are given by 

^:w=<(/i)+4^ia(/i)c,w, 



The anomalous dimension tensor 7^ governs the mixing of the double insertion of Q± 
and Qa into Q u (see Fig. EJ), while 7± describes the self- mixing of Q± and was computed 
in [29]. The anomalous dimensions read: 



72? = 2(1 ±3), 



P(e) 

7±V 



0. 



P(es) 

7± * 



52 



1 ±3). 



We have defined the matrix 7^ as 



p(k) 

i±y 



2 7±,„ 



sin 6> w 7 ± 



V(k) 



(31) 



(32) 



with the superscripts A and V denoting the contributions stemming from double insertion 
of (Q±,Q 9 A ) and (Q±,Q V ), respectively. The LO result agrees with [5,14]. The other 
contributions are new. 
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Figure 6: Sample diagrams for the NLO mixing of Qa and Q± into Q h 



The anomalous dimension of Qa vanishes and the RGE in the penguin sector is the 
linear equation 



where 



'4C + {p)C A 
Cp(fi) = I 4C-(ji)C A 



T 

Ip 




(33) 



(34) 



,0 7 „ 

The RGE for single insertions can be solved explicitly using the method described in [30]. 



2.3 Below \x c 

At fi c , i.e. the scale of the charm quark mass, the charm quark is integrated out and 
removed as a degree of freedom. All necessary matrix elements are given in [5,14] - no new 
contributions arise to the orders considered here. There are some new terms stemming 
from the expansion of m c (/z c ) about m c (m c ) in these expressions, though, and we collect 
these results for convenience. 

The matching in the box sector leads to the following matrix elements: 



rf ( Mc ) = 5 + T -^- In x T + 4 In 



1 — x T 



mi 



(35) 



where x T = m T jm c and m c = m c (/i c ). Neglecting the lepton masses for the electron and 
muon, the above formula yields 



rfPiHc) = 5 + 41n 



ml 



(36) 



We have defined the matrix elements for lepton flavour I by 

(Of 0O> = ^r?UM(Q y )® 

where (Qf (/i c )) denotes the double insertion of the operators in the box sector. In the 
penguin sector we find 



(37) 



r£ (1 V c ) = (l±3)(l-ln 



mi 



(3* 



where 

m^)) = ^r P ^\^)(Q^\ (39) 
and (Q±([/, c )) denotes the double insertion of the operators in the penguin sector. 

2.4 Final Analytic Expression for P C (X) 

Now all that remains to do is to combine all relevant terms and compute the box and 
penguin contributions to the function X l (x c ) defined in Eq. (Tj[|). Here we closely follow [14]. 
Let us start with the box contribution. We expand the result as 

c ^ - ^ (^°>m + w + w) (40) 

and express the running charm quark mass m c (/x c ) in terms of initial condition m c (m c ), 

„ w (i + a -^e + + £<->) *«k), (4i) 

/ (0) ,a \ 

where we defined k c = rjc and ?7 C = a s (fi c )/a s (m c ), with the individual contributions 

d e) = ^c-i), 

Po 



-,(") -.(0) o ..(e) 



/m / Jm p\ Jin \ , \ //r,\ 

+ ^(ir~7%J (1 -"< (42) 

We find the following expansion coefficients for C l B : 
Cf\» c ) = C^ c ), 

Cf\» c ) = C^ c ) + CWQiltf + 4CSV c ) 2 pf (e) (/ic), 
+ 4^V c )Vf (es Vc) + 4Cj?0O V (e) (/. c )d 1} 

+ SC^i^C^i^pf (1) (// c ) + 4^ ) ( / u c ) 2 pf (1) (Pc)d e) - (43) 



10 



We obtain the parameters pf by inserting the expansion of m c (/i c ) into the expressions for 



rf (see Sec.lOl: 



X -y 



In 



X T [1 K c 

1 — X T 



lnx r 



pf ^ = 



n B{es) 
Ft 



4K C e 



k c — x T ( 1 — In ^ L 



(44) 



The corresponding expressions for the electron and the muon, where we can neglect the 
masses, are given by 



(45) 



The penguin contribution to the function X l (x c ) can be obtained in the same way. 
Expanding the Wilson coefficients Cp(p c ) as 



Cp(fl c ) = K c 



( 4-C? W + + -^CPW) . (46) 



32 \a s (fi c ) 
we find the following contributions: 

4V) = ^°>(/z c ), 

i=± 

c ( p\p c ) = c^ s \p c ) + c^\p c )e ] + a p(i) (/i C )d e) + a p(0 Vc)d es) 

+ 4£ U ies \p c ) + pf (e) (/, c )d 1) + P^\p c )6 e) ) cjpMcfM 



i=± 



+ 4$>f (1) (p c ) (ct\^)Cf{» c ) +Cf (p c )4 e) (^c) 
i=± 

+ 4^pf (e) (p c )Cf ) (p c )4° ) (^). 



(47) 



i=± 



Again we obtain the parameters pf by inserting the expansion of m c (p c ) into the expres- 
sions for rf : 



p p(i) = r ^a)( mc ) + (i±3)ln Kc 
P? e) = o, 

p p j es) = (i±s)e c e) . 



(48) 
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M w 


(80.403 ± 0.029) GeV 


[31] 


ot s (M z ) 


0.1176 ±0.0020 


[31] 


M z 


(91.1876 ± 0.0021) GeV 


[31] 


a(M z ) 


1/127.9 


[31] 


M t 


(172.6 ± 1.4) GeV 


[32] 


sin 2 ^ s 


0.23122 ±0.00015 


[31] 


m b {m b ) 


(4.164 ±0.025) GeV 


[33] 


Gp 


1.166 37 x 10- 5 GeV"^ 


[31] 


m c {m c ) 


(1.286 ±0.013) GeV 


[33] 


X 


0.2255 ±0.0007 


[34] 


M H 


(155 ±40) GeV 




\Vcb\ 


(4.15 ±0.09) x 10~ 2 


[35] 


m T 


(1776.99l£aj) MeV 


[31] 


p 


U.l / ±l_o .017 


[35] 








fj 


0.343 ±0.016 


[35] 



Table 1: Input parameters used in our numerical analysis. 



The final result for X 1 is then 

X l (x c ) = C P (nc) + C l B ( f i c ). (49) 

The corresponding expressions for Cp(/z c ) and C l B {ji c ) can be found in Eqs. (140!) and (1461) . 
respectively. Eq. §5§ then yields the contribution to the branching fraction. 



3 Final Results and Numerical Discussion 

Having all necessary ingredients at hand we will discuss the numerical implications of our 
results, where we use the input parameters given in Tab. HJ Our numerical procedure 
follows closely the one of Ref. [14]. In particular we use the numerical solution of the 
RGE of the program RunDec [36] to compute a s (/i c ) from a s {Mz) and neglect all terms 
proportional to (3 es . We have checked numerically that this is indeed justifiecll. 

The dependence of P C {X) on the parameter fi c can be seen in Fig. [7J We use central 
values for all relevant input parameters of Tab. [T] and fix /i& = 5 GeV and = 80 GeV. 
The dashed line shows P C {X) as a function of /x c including the NNLO QCD corrections, as 
computed in [14] where the parameter x c equals m 2 /M^. The dashed-dotted line shows 
the same quantity, but using our improved definition of x c , see Eq. ([3]). We observe that 
this line is shifted by about 0.5% compared to P C {X) using the conventional definition 
of x c . The dotted and the solid lines show the results including LO QED and the NLO 
electroweak corrections, respectively. We see that including the full electroweak corrections, 
P C {X) is increased by another 1.5% as compared to the pure NNLO QCD result with the 
improved definition of x c . Also the cancellation of the scheme dependence between the LO 
QED and the NLO electroweak contribution is clearly visible. 

The explicit analytic expression for P C {X) including the complete NNLO corrections 
is so complicated and long that we derive an approximate formula. Setting A = 0.2255 
and mt(mt) = 163.0GeV we derive an approximate formula for P C {X) that summarises the 
dominant parametric and theoretical uncertainties due to m c (m c ), a s (Mz), /i c , p>w, and 

2 We thank Ulrich Haisch for providing us with his program for the QED running of a s . 



12 



0.39 



0.38 - 



OT 0.37 - 



0.36 - 



0.35 




_L 



_L 



_L 



1.25 1.5 



1.75 2 2.25 2.5 2.75 3 

Mc [GeV] 

Figure 7: P C (X) as a function of fi c at NNLO QCD (dashed dotted line), including LO QED 
(dotted line), and NLO electroweak corrections (solid line). The dashed line shows P C (X) at 
NNLO QCD where the definition x c = m c /M]y is used. 



fib- It reads 



P C (X) = 0.38049 



m c {m, 



0.5081 



1.30GeV/ 



± 0.008707 



mJm C/ 



1.30GeV7 



0.5276 



where 



In 



m c (m c 



1.30GeVy 



In 



oc s {M z ) 
0.1176 

a s (M z ) 
0.1176 



QtsjMz) 
0.1176 



1.0192 



1.8970 



K- T l J 3 

""tJ^tncOla 



h3 



(50) 



LI 



1,3 



(51) 



and the sum includes the expansion coefficients and given in Tab. [3 The above 
formula approximates the central value of the full NNLO QCD result plus electroweak 
corrections with an accuracy of ±0.05% in the ranges 1.15 GeV < m c (m c ) < 1.45 GeV, 
0.114 < a s (M z ) < 0.122, while the scale uncertainty for varying 1.0 GeV < fi c < 3.0 GeV, 
40 GeV < n w < 160 GeV, and 2.5 GeV < /i b < 10.0 GeV is correct up to ±2.3% in 
Eq. ( ISTfl) . The uncertainties due to m t (m t ), and the different methods of computing a s ([i c ) 
from a s (Mz), which are not quantified above, are all below ±0.2%. For A = 0.2255 we 
find P C (X) = 0.372 ±0.015, where 42% of the error are related to the remaining theoretical 
uncertainty and 58% to the uncertainties in m c and a s . In the future one could utilise the 
correlation of m c and a s in Ref. [33] to further reduce the parametric uncertainty. 
Finally we provide an updated number for the branching ratio: 



B (K + -> 7T + z/z/(7)) = (8.51 ± 0.20 ± 0.36) x 10 



i-ii 



(52) 
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kio = 1.6624 


kqi = —2.3537 


k u = -1.5862 


K2Q = 1.5036 


k 02 = -4.3477 


eio = -0.3537 


e 01 = 0.6003 


en = -4.7652 


e 20 = 1.0253 


e 02 = 0.8866 



Table 2: The coefficients Kij and arising in the approximate formula for P C (X). 



The first error stems from the uncertainties in the CKM parameters. The second error is 
related to the uncertainties in m c , m t , and a s , where all three quantities contribute in equal 
shares. The dependence on Mh is completely negligible (below one per mil). The last error 
quantifies the remaining theoretical uncertainty. Here the main contributions stem from 
the uncertainty in <5P CjM and X t , where we used an error of 2%. In detail, the contributions 
to the theory error are («+ : 6%, X t : 38%, P c : 17%, 5P CtU : 39%), respectively. All errors 
have been added in quadrature. 

4 Conclusion 

In this paper we have calculated the 0(a) and 0(aa s ) anomalous dimensions and the 
electroweak matching corrections of the charm quark contribution relevant for the rare 
decay K + — > n + uu. The parametric dependence of the relevant parameter P C (X) plus its 
theoretical uncertainty is summarised in an approximate but very accurate formula. 

P C (X) is increased by up to 2% as compared to the previously known results [14]. This 
change is of the same order of magnitude as the remaining scale uncertainties after the 
NNLO QCD calculation. Together with the recently achieved very precise determination of 
the hadronic matrix elements [17], further improvements on the long-distance contribution 
of the charm quark [20], and the complete electroweak matching corrections for the top 
quark contribution [37] the theoretical prediction of the branching ratio B(K + — > 7i + uu) 
will reach an exceptional degree of precision, with the uncertainties mainly due to the 
CKM parameters. 

The latter errors will be reduced in the coming years by the 5-physics experiments and 
a precise measurement of the branching ratio B(K + — > 7T + uu) will provide a unique test 
of the flavour sector of the SM and its extensions. 
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